Lowest Landau Level vortex structure of a Bose-Einstein condensate rotating in a 

harmonic plus quartic trap 
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We investigate the vortex patterns appearing in a two-dimensional annular Bose-Einstein con- 
densate rotating in a quadratic plus quartic confining potential. We show that in the limit of small 
anharmonicity the Gross-Pitaevskii energy can be minimized amongst the Lowest Landau Level 
wave functions and use this particular form to get theoretical results in the spirit of [A. Aftalion 
X. Blanc F. Nier, Phys. Rev. A 73, 011601 (R) (2006)]. In particular, we show that the vortex 
pattern is infinite but not uniform. We also compute numerically the complete vortex structure: it 
is an Abrikosov lattice strongly distorted near the edges of the condensate with multiply quantized 
vortices appearing at the center of the trap. 

PACS numbers: 03.75.Lm,67.85.Bc 
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I. INTRODUCTION 

When rotated, superfluids are expected to exhibit a 
rich structure of quantized vortex states. This has been 
experimentally confirmed for Bose-Einstein condensates: 
the techniques developed to create a single vortex state 
[U, rapidly led to arrays containing up to several hun- 
dreds of vortices 0, 13, S 111 ■ Typically these condensates 
are confined by harmonic potentials, and the large vortex 
arrays are obtained for angular velocities f2 approaching 
the radial trap oscillator frequency uj. 
With these type of potentials, the limit — > w is singu- 
lar: as the centrifugal force balances the trapping force 
the Thomas-Fermi radius of the condensate diverges and 
the central density decreases towards zero. 
In the experiments 0, [3 , a blue detuned laser directed 
along the axial direction adds a quartic component to the 
usual harmonic potential, resulting in a potential of the 
form : 



V(r) = -niLj'^r'^ + -r^ 
2 4 



(LI) 



where m is the mass of the atoms. The nice feature of 
this potential is that the centrifugal force, varying as 
Q?r can always be compensated by the trapping force 
varying as — (mw^r -t- kr^), so that one can explore the 
region 17 > w. 

Condensates rotating in this kind of trap present a 
very rich variety of vortex phases. At sufficiently slow 
rotation speeds, the system is not expected to have any 
vortex. For increasing SI there is a sequence of states 
with more and more vorticity, where singly and multiply 
uantized vortices may be present at the same time 

"Bp. 

At rotation speeds ~ (jj, a large triangular array 
of singly quantized vortices appears, which is not 
qualitatively different from what happens in a purely 
harmonic trap. The difference is when U, is taken larger 
than uj: the condensate continues to expand radially 
and there is a critical speed fie at which the centrifugal 



force creates a hole in the condensate at the center of 
the trap [H, [13, [3 ■ For rotation speeds > fie there is 
a vortex lattice in an annular region and a central hole 
where the density vanishes and around which there is a 
circulation. However, this hole is not usually referred to 
as a giant vortex as not all the circulation is contained 
in a central vortex (T^ . 

Finally, for even larger rotation speeds the condensate 
continues to expand radially with constant area so that 
the width of the annulus decreases and all the vortices 
are expected to retreat into the central hole, resulting 
in a pure irrotational state with macroscopic circulation 
[Ulillill, often called a giant vortex [illI3|. 
In this paper we address the regime > fie- We would 
like to know the precise repartition of vortices, both 
in the annular Thomas-Fermi region where we expect 
that they form a regular lattice, and in the region of 
low density. Indeed, we believe that there should be 
"invisible" vortices beyond the external radius of the 
condensate (as is the case for a harmonically trapped 
condensate in fast rotation and in the central hole 

In the limit of weak anharmonicity, the single particle 
states are restricted to the Lowest Landau Level (LLL) 
wave functions fl3|. Thus our analysis consists in mini- 
mizing the Gross-Pitaevskii energy in the LLL. Using the 
explicit expression of the projection onto the LLL we get 
theoretical results in the spirit of [l^, [l^ . We prove that 
a wave function minimizing the Gross-Pitaevskii energy 
in the LLL cannot have a finite number of zeroes, thus 
the vortices of the condensate cannot lie in a bounded 
domain. This implies that there should be an infinite 
number of invisible vortices. We also construct critical 
points for the Gross-Pitaevskii energy by distorting a 
regular hexagonal vortex lattice in the following way: 
first, using a theta function as is done for the Abrikosov 
problem ^20i] . we construct a function u which modulus 
is periodic over an infinite regular hexagonal lattice and 
which vanishes only at the points of the lattice with 
simple zeroes. Next we multiply w by a slow varying 
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profile a (which corresponds to the square root of the 
coarse^rain average of the atom density over several 
cells |14l. l2l| \ . We then project au onto the LLL. It can 
be shown [19] that projecting this function onto the LLL 
implies a distortion of the vortex lattice. 
The state of a condensate in the LLL is entirely known 
from the positions of its vortices. This allows us to 
adapt the numerical method of and minimize the 
Gross-Pitaevskii energy functional of the condensate by 
varying the locations of its vortices. As is expected in 
the range of parameters we chose, the condensate lies 
in an annular domain in which vortices are regularly 
distributed j^, il3|] . We also find that there are invisible 
vortices both beyond the outer radius of the condensate 
and in the central hole. The vortex lattice is strongly 
distorted at the inner and outer edges of the condensate 
and we find that multiply quantized vortices appear at 
the center of the trap. 

The paper is organized as follows: in Section we 
describe our model and general approach. We show 
how reducing the minimization of the Gross-Pitaevskii 
energy to the LLL leads to a new variational problem 
for which we derive an Euler-Lagrange equation satisfied 
by the minimizer. We then use this equation in Section 
lllll to show that the minimizer necessarily has an infinite 
number of zeroes and to construct critical points for 
the energy. In Section IIVI we describe our numerical 
method and comment our results, relating them to our 
theoretical approach. Finally we give our conclusions in 
Section |Vl 



II. MODEL-APPROACH 

A. Gross-Pitaevskii energy and LLL reduction 

We assume a strong confinement along the rotation 
axis, so that we are in an almost 2-D situation. We take 
cj, huj and ^tijiraiS) as units of frequency, energy and 
length respectively and consider the Gross-Pitaevskii en- 
ergy of the condensate in the rotating frame. 



^GP(^) 



l-r!2 



dx, (II.l) 



where x = {xi,X2), = {—X2,xi) and G is a dimen- 
sionless coefficient characterizing the strength of atomic 
interactions (see [13, [13] for example). We minimize £gp 
under the mass constraint / IV'P = 1- The main idea of 
our analysis is to restrict the minimization of £gp to the 
first eigenspace of the operator — (V — iflx-^) , corre- 
sponding to the eigenvalue 11. This is the Lowest Landau 
Level, introduced in the context of Bose-Einstein con- 
densation by Ho [2l[ and then successfully used in differ- 



ent studies on harmonically trapped condensates rotat- 
ing at speeds close to the trap frequency (see for example 



U, 122, 123(). In [i, [13] the LLL approximation is used 
to determine the phase diagram of a condensate trapped 
by a harmonic plus quartic potential in the limit of weak 
interactions and small anharmonicity. 
The LLL consists of functions of the form ip{z) = 
/(z)e~^l^l where z is the complex variable xi + 1x2 
and / is a holomorphic function. For an LLL function -0, 
the Gross-Pitaevskii energy reduces to 



£lll{'>P) = 



G. 



2 |zp + -|z|M|^|2 + _|^|4jrf^^ 



(II.2) 



so that we will minimize the energy £lll in the LLL 
under the mass constraint / JV']^ = 1- 
Minimizing pi.2p without the assumption that ip is in 
the LLL gives the main scales of the problem. One gets 
the Thomas-Fermi distribution 



\^TF{r)\' - max ^+ ^ J 



(II.3) 



where r is the radial coordinate. We refer to [12, [3] 
for the detailed study of such a distribution. In particu- 
lar one gets the expression of fie (see equation piL19|) ). 
which allows to deduce that for ]1 — fi] at most of the 
order of k'^/^G^/^ we have O ~ flc- Then £lllW has 
the order of k^^^G^/^ and the spatial extension of V'tf is 
of the order of {G/kf^. 

The LLL approximation is valid when two conditions are 
fulfilled [lj,|2l|: («) the excess energy £lll{iPtf) is small 
compared to the splitting between the LLL and the first 
excited Landau level (~ 2 in our units), and (m) the 
spatial extension of the condensate is large compared to 
the vortices interdistance, so that coarse-grain averages 
of atomic and vortex densities over several vortex cells 
make sense. Using the scalings derived from (jll.3|) . this 
conditions reduce to 

fc<G< ^ and ll-fl] <fc2/3Gi/3. (II.4) 
Vk 



In Section IIII CI we will give more details on how the 
scalings pL4p are derived and see that taking into ac- 
count the LLL constraint does not modify the scales of 
the problem. Indeed, we show in Section IIII Bl that al- 
though a function such as pi.Sp cannot be in the LLL, a 
proper distribution of vortices allows one to approximate 
very closely such a profile in the LLL, modifying only 
slightly the averaged atomic density of the condensate. 



B. Mathematical framework for the small 
anharmonicity regime 

We now describe the mathematical framework we are 
joing to use in Section IIIII It was introduced in [l^. 
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|19| for the theoretical study of an harmonically trapped 
condensate. 

We define a small parameter 

corresponding to our small anharmonicity regime and 
study the asymptotics of the problem as e — > 0. We 
take 

n = l + f30^ (11.6) 

with (3 ~ G^^^ so that we are in the range of pa- 
rameters (jll.4[) where the LLL approximation is justi- 
fied. We rescale distances by making the change of vari- 
ables (j){z) = ^/eip{y/ez), the energy becomes Slll{iP) = 
£Elll{4>) with 

ELLdcp) = ^ [[~P\z? + W + ^101^) dz. 

(n.7) 

We minimize this new energy amongst all functions (p in 
the LLL satisfying the mass constraint / |0p = 1. 
For every (p in the LLL we define the function f{z) = 
(/)(z)el^l By definition of the LLL, / belongs to the 
Fock-Bargmann space 

J^e = |/ holomorphic |/pe-l^l'/^dz < c»| . (ILS) 

The space J> is a Hilbert space for the scalar product 

{f,g)= I 7M.g(z)e-l^l'/^dz. 
Jc 

The point of introducing such a space is that the orthog- 
onal projection of any function g onto J-^ is explicitly 
known [M,15]: 

Tl,{g){z) = — / e^-"/^e-l^'l'/^5(^')rf^', (n.9) 
Jc 

so that if we write the variational problem (|TL7} for / 
rather than for 0, namely if we look for / G J-^ minimiz- 
ing 

+ (11.10) 

under the constraint (/. /) = 1, we are able to derive that 
/ satisfies the equation 

- PMJ + \ {M^f + sMj) + Gne(e-I-I'/^|/|2/) = f^f 

(11.11) 

where fi is the chemical potential coming from the mass 
constraint and Mg is the operator defined by 

M, = ed,z. 



Indeed, the weak derivative of the energy Flll at / along 
g is given by 

+f l/P/ge-^i^i^/^) dz. 

Using integration by parts and dgf ~ d^g = on the first 
term and the fact that II^ (g) ~ g on the second term, we 
obtain pi.lip . 

We give two equivalent forms of (jlLlip that we shall need 
in the sequel: 

- me! + \ (Me / + £M,/) + ^/(ea,)[/2(z/2)] = 

(11.12) 

(-/3 + l£)n,(|zp/) + |n,(|zpn,(|zp/)) 

-f Gn,(-i^i'/^n/p/) = M/. (11.13) 

The operator f{edz) is defined by 

if f{z) = ^Cfcz'^. Equation pi.l2p is obtained 
from (jlLllj) as in [l^ with some algebra on the 
non-linear term. To get pi.lSp we use an integra- 
tion by parts to show that M^/ = Hgdzp/), then 
M^f = Bedzpnedzp/)) and we get the result. 
In the following section we prove our theoretical results. 
Using pi.l2p we show below that any minimizer / 
of pi.lOp (and thus any minimizer (j> of (HLZi) has 
an infinite number of zeroes for s small enough. The 
construction of critical points is done using pi.lSp . 



III. ANALYTICAL STUDY 

In this section we present the results we are able to 
derive from equations pi.l2p and pi.lSp . 

A. Infinite number of zeroes 

We show that any minimizer / of pi.lOp (and there- 
fore any minimizer ip of pi.2p ) has an infinite number of 
zeroes if e is small enough. The argument is by contra- 
diction and in two steps. 

1. Suppose / has a finite number of zeroes. Then one 
may write /(z) = F(z)e'^(^) where P is a polyno- 
mial and is a holomorphic function. Now f € J^g 
and the condition |/pe~l^l ^^dz < oo implies 
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that Re{(f{z)) < \z\'^/{2e). It is well-known (see 
[26| for example) that a holomorphic function can 
satisfy this condition only if it is a polynomial of 
degree less than 2. Therefore we know that 



/(z) = P{z)e 



(III.l) 



and the integrability condition on / implies a2 < 
l/(2e). Injecting (fllLTTl in ([IlTT|1 and comparing 
the exponential growth of the different terms of 
(|IlTT|) as in ^ yields ai = as = 0. So, if / 
has a finite number of zeroes, it is a polynomial. 

2. Now, suppose / is a polynomial of degree 
n and inject this in (jII.12p . The term 
-/3M,,/ + i (M|/ + eAf. /) is a polynomial of de- 
gree n, therefore (jII.12p implies that the term 
Y/(e9z)[/2(z/2)] is also of degree n. But 
[edz)^\P{z is of degree 2n — fc, so that / must 
be of the form /(z) = cz". Injecting this a last time 
in (|II.12p . using the improved Stirling [l^] formula 
and the condition (/, /) = 1 yields a condition on 
n: 



lJ.-\- Pe> —fine - 



27reY^ 



(III.2) 



We now bound the chemical potential /i: taking 
the JFe-scalar product of each side of (|II.lip with / 
yields 

tJi<2FLLL{f) + fi^~^+'^. (III.3) 

Here we used the fact that the spectrum of —(3M^ + 
J (Mg -I- eMg) is bounded below uniformly with 
respect to e. We shall see in Section IIII CI (see 
(|III.22p i that 



FlllU) = ELLLife 



5 V stt y 



Thus we have 



If! 

16 



> -l3ne + 



2Trey/n 



(III.4) 



and the left-hand side of pil.4p is bounded uni- 
formly with respect to s. Minimizing the right- 
hand side of pil.4p with respect to n (taken as a 
continuous variable as it should be very large when 
e is small) for fixed e yields 



(III.5) 



and 



6 fibCy^^ _ ^2 _ ^ _ 



5 V Stt y 2 16 

5 /Ge'i/i2x4/5 



> - , 

- 4 V 27r 



-1/5N 



(IIL6) 



which is a contradiction if e is small enough. 

In the rest of the paper we take an e small enough for the 
inequation (|III.6|1 not to be verified. We then conclude 
that / has an infinite number of zeroes. As there is a limit 
on how close vortices can be in this regime (see [HI), they 
cannot lie in a bounded domain and the vortex pattern 
extends to infinity. 



B. Construction of critical points 

As is the case for harmonically trapped condensates 
in fast rotation, we expect that in the range of parame- 
ters we explore the scales of the problem will decouple. 
Namely we expect any minimizer of (jll.7[) to be of the 
form au where u varies on the scale of the vortex pattern 
(which is small compared to the size of the condensate) 
and a is a slow varying profile giving the general shape 
of the condensate. In this section we show that although 
such a function is not in the LLL, one can approach it by 
an LLL function which is an almost critical point for the 
energy (jll.zp . We also relate the effect of the projection 
onto the LLL to the distortion of a regular lattice in the 
region of low atomic density. 
More precisely we introduce 

Uriz) = e-l^l'/2^/,(z), Mz) = e^V^-e (^^^z, 

(111.7) 

where t = t/? + irj is any complex number and 



e(«. 



+00 

7 Y. 



(_l)"e"'^("+l/2)^g(2n+l)7rTO^ (III.J 



The Q function has the property (see |28j for more de- 
tails) e{v + k + It,t) = (-l)'=+'e-2*'^'"e-"'^e(w,r) so 

that |'iir(z)| is periodic over the lattice y^^Z© ^J^Zt, 

and vanishes at each point of the lattice. 
The interest of introducing such functions is twofold. 
Firstly it is known (28| that any function v whose modu- 
lus is periodic over the lattice \J^'^ © \J^'^'''^ vanishes 
exactly on the points of the lattice with simple zeroes 
and such that g = we'^' is holomorphic must be pro- 
portional to Ur- Secondly, the function /t is a solution 
to the Abrikosov problem (see [H, [H [lo, Hi] ) 



n(IM 



It) = Kfr, with = (lUrP) &(t). 



(IIL9) 
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and 



&(r) 



Iri 



(III.IO) 



kjei, 



Equation (|III.9[) is similar to pi.l3|) without the potential 
term and with fi = so that one can expect to obtain 
a solution of pi.l3p by a slight modification of fr- We 
refer to [l8[ and the references therein for details on the 
quantity 6(t). Let us just mention that it is minimum 
(6(t) ^ 1.16) for T = e^''^/^, which corresponds to a 
hexagonal lattice. 
We define 



(n,(a/r),n,(a/,)) 



1/2 



and 



(III.ll) 



(III.12) 



where a € C°'^/^(C, C) is a slow varying profile with 
compact support and J |ap — 1. We use Lemma 5.4 in 
[191, which states that 



Ua.r — CtUr + 0{e 



1/4N 



(III.13) 



This, together with the fact that |wrP is periodic over a 
lattice of period bounded by 0(e^/^) yields 

(III.14) 



and 



FLLLifa.r) — £'lLl(Wq,t) — 



(III.15) 



Equations piI.14p and piI.15|) are the scale decoupling 
we expected: the only contribution of the lattice to the 
energy is through the coefficient 6(t), which is minimum 
for the hexagonal lattice, r = e^'^'^l'^ . From now on wc 
shall note h — 6(6^"/"^) ~ 1.16. We now have to minimize 
pil.151) . which is an energy similar to pi.2p . with respect 
to the profile a which is not in the LLL, so that the 
minimization is performed in the usual way. We obtain 
the Thomas-Fermi profile 



\a\\z) 



bG 



,0 



(III.16) 



where v is the chemical potential associated with the con- 
straint /lap = 1 and b takes into account the vortices 



contribution. With such a profile and taking fi = v in 
piI.14p we get 

+ l£)n,(|z| + Jn, (izpnakP/.,.)) 

+ Gn,(e-l^l'/^|/„,,|2/a,r) = JyU,r + 0(£^/") (III.17) 

so that fa^T is almost a critical point for the energy 
(IlLTOl) . 

We now explain the effect of our procedure on the vortex 
pattern. The property Ua^r = ait^r + 0(e^/^) shows that 
on the support of a the zeroes of distributed 
close to those of Ur, on a regular hexagonal lattice. To 
get information on the "invisible vortices" lying in the 
region of low density, we evaluate the number N{R) of 
zeroes of fa^r in a ball of radius R. The Cauchy formula 
yields 



R 



2Tr 



N(R) = d9 - 

27re7o /e«^'/ee-kT/ea(z')/r(2'e-*«)dz' 

(III.18) 

Using a Laplace method to evaluate the integrals, we 
see that their ratio is bounded for large R, so that 
N{R) oc R/e. A regular lattice would give N{R) cx R'^/e, 
so we deduce that the lattice is strongly distorted out- 
side the external radius of the condensate. Note that this 
method does not give information on the distribution of 
vortices in the central hole when the condensate has an 
annular form, for this would correspond to R small, and 
the Laplace method is not efficient in this case. We refer 
to our numerical simulations in Sect ion HVl for the vortex 
structure in the central hole. 



C. Evaluation of physical quantities 



Here we analyze the TF profile pil.161) . or rather 
ct{z) = y/e ^ a[z^ ^), to evaluate some relevant physi- 
cal quantities in the original scaling. 
The analysis of a profile such as piI.16p has already been 
done in [l3 . [3] , so we only adapt and summarize the re- 
sults. The critical rotation speed for the condensate to 
develop a central hole is 



\ 8tt 



1/3 



(III.19) 



For subcritical velocities, the behavior of the condensate 
is not qualitatively different from that of a harmonically 
trapped condensate, so we focus on velocities il > ^Ic 

which is equivalent to (3 > Then the inner 

and outer radius of the condensate R± are given by the 
relations 



R% + Rt^ 4(3k 



-1/3 



Ri - Rt 



1/3 



(III.20) 
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the chemical potential is 



trapped condensate) would be to take 



(III.21) 



and the energy is 

This is slightly different from the result obtained in [l^l 
with the solid-body approximation. Indeed, our analysis 
allows one, through the coefficient b, to take into account 
more precisely the contribution of the vortex patterns to 
the energy and atomic density, whereas the solid-body 
approximation implies an infinite regular lattice with pre- 
scribed volume of the cell. These results and the fact that 
& ~ 1.16 imply that our analysis is valid within the range 
of parameters we announced in Section 2.1 : equations 
pil.20p yield that the spatial extension of the condensate 
is large compared to 1 (which is the order of the vortex 
pattern spacing) if 

We now see that the condition f3 ~ G^/^ (or \1 - n\ < 
^2/3(^i/3-j jg equivalent to fl ^ Qc and ensures that the 
energy Slll is of order fc^/'^G^/'^, so that the energy con- 
dition for the LLL approximation to be valid reduces to 



G< 



Now, as we have both k <^ G and G ^ we see that 
necessarily <C 1, which justifies our study of a small 
anharmonicity regime. 



IV. NUMERICAL SIMULATIONS 

A. Approach 

We want to numerically approach the minimizer of 
Slll in the LLL. We write 



(z) = P(z)e-"l^l'/2 



(IV.l) 



with P a holomorphic function. As polynomials are dense 
in it is reasonable to fix an integer n and to restrict 
the analysis to functions 4> where P is a polynomial of 
degree less than n. We write our trial functions as in 
p^ . where A = \\(l)\\j^2^'^ is the normalization factor: 

n 

(^(z) = A[|(z-Zj)e-"l"l'/2 (IV.2) 

and vary the locations Zj of vortices. An alternative 
method (used for example in [22] for a harmonically 



J „-n\z\^/2 



and vary the coefficients bj . The interest of our approach 
is to give a direct access to the exact repartition of vor- 
tices, whereas the alternative method would require to 
compute the roots of a polynomial of degree n, which is 
a delicate task for large n. In particular, varying the co- 
efficients could probably not give the precise locations of 
invisible vortices. 

We used a conjugate gradient method with a Goldstein 
and Price line-search. The integrals are computed using 
the Gauss-Hermite method, and we take enough Gauss 
points for the computations to be exact. This results 
in quite expensive calculations, but we have been able 
to numerically construct condensates with up to ~ 120 
vortices. 



B. Results 




FIG. 1: Vortex structure and atomic density for G = 3, fc = 
10"", /3 = 1 (n - 1 = 2.2 10"^). There are 67 vortices in total, 
the central vortex is constituted of 11 single vortices. 

We show in Fig. [T] and Fig. [2] typical examples of 
configurations we numerically computed. The qualita- 
tive features of the vortex patterns and atomic densities 
confirm our theoretical results and are in good agreement 
with existing theoretical and numerical studies [H, [3] ■ 
Note however that the numerics become quite intricate 
for large number of vortices, which accounts for the rel- 
ative lack of symmetry of Fig. [21 

As was expected, the condensate develops a central 
hole and visible vortices are regularly distributed in the 
annular region of significant atomic density. Our compu- 
tations also show a distortion of the vortex pattern near 
the external radius of the condensate, resulting in invis- 
ible vortices in the exterior region of low atomic de nsity 
as is the case for harmonically trapped condensates [ij] . 
Some vortices also lie in the central hole as theoretically 
predicted (see for example jT^ ) and we can get informa- 
tion on their precise locations: we observe a distortion of 
the regular lattice near the inner radius of the conden- 
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FIG. 2: Vortex structure and atomic density for G = 3,k — 
10"^/3 = 1 (fi - 1 = 4.610"''). There are 119 vortices in 
total, the central vortex is constituted of 20 single vortices. 




cantly. We find good agreement of our numerical results 
and analytical study, with energies typically differing by 
~ 10"^ to ~ 10"^ in relative value. 



52 54 53 5S 
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FIG. 4: Minimum energy as a function of the number of vor- 
tices in the trial wave function (G = 3,fc = 10"^, (3 = 1). 



CONCLUSION 



FIG. 3: Two example of vortex configurations for G — 3,k — 
10""*, /3 = 1 - 1 = 2.2 10"^), respectively with n = 60 and 
n = 67 vortices. 

sate, resulting in isolated singly-quantized vortices encir- 
cling a central multiply quantized vortex. We computed 
configurations for which this central vortex has up to 20 
units of circulation while there is a total of 32 units of 
circulation in the entire hole and 83 visible vortices. The 
number of vortices (both visible and invisible) increases 
with increasing 17 or /3, but since our scaling does not 
allow to explore a large domain of fl when G is fixed, 
we mainly varied k. The total vorticity of the system 
increases with decreasing k. 

All vortices do not have the same contribution to the 
energy: as n increases, the vortex pattern in the annu- 
lar region of significant atomic density remains the same 
up to possible rotations, with the additional vortices first 
gathering in the central vortex, then constituting the dis- 
torted lattice near the inner boundary of the condensate 
and finally occupying the distorted sites beyond the ex- 
ternal radius. With increasing n (see Fig. SJ, the energy 
reaches a first plateau when the central figure is formed, 
constituted of the visible vortices and the vortices in the 
central hole. A second plateau is reached when enough 
distorted sites beyond the external radius are occupied. 
For example, for the parameters corresponding to Fig. [1] 
and Fig. [3](fc = 10^**, G = 3, (3 = 1), the energy varies by 
~ ±10~^ in relative value when n is increased from 60 to 
67 (Fig. [3l) and the atomic density does not vary signifi- 



We have investigated the vortex structure of a two di- 
mensional annular Bose-Einstcin condensate rotating in 
a quadratic plus quartic confining potential. We focused 
on the regime where the state of reference is a vortex 
lattice encircling a central hole carrying a macroscopic 
circulation. We developed a theoretical and numerical 
approach based on a small anharmonicity regime and a 
reduction to the Lowest Landau Level states. Using a 
theoretical method developed in [l^, [l^ we showed that 
there is an infinite number of vortices in the condensate 
so that they cannot lie in a bounded domain. We then 
analytically and numerically constructed critical points 
for the reduced Gross-Pitaevskii energy and get further 
information on their vortices. We find that they are reg- 
ularly distributed on a hexagonal lattice in the annular 
Thomas- Fermi region where the atomic density takes sig- 
nificant values, and that this pattern is strongly distorted 
at both edges of the annulus. In particular we find nu- 
merical evidence of multiply quantized vortices appear- 
ing at the center of the trap. Our results agree with 
and complete existing studies [12, [13, [III , showing how 
the distortion of the vortex pattern modifies the results 
obtained with the solid-body approximation. 
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